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A  me cha r> i 1  analog  Vertical  Energy  Absorption  Model  (VRAM)  is  developed  to 
predict  the  dynamics  of  an  Air  Cushion  Landing  System  (ACL3)  in  ebe  vertical 
dimension.  Three  degrees  of  freedom  and  thus  three  primary  modes  of  oscillation 
arc  investigated i  heave,  pitch  and  roll.  Data  from  Air  Force  Flight  Dynamics 
Laboratory  Tests  of  a  full-scale  Australian  Jindtvik  drone  are  used  to  develop  and 
verify  the  model. 

The  VEAM  study  demonstrates  that  the  use  of  a  mechanical  analog  prediction 
scheme  for  an  Air-Cushion  Landing  System  has  sufficient  merit  to  warrant  further 
invest ige.t:  on.  It  shows  that  a  mechanical  analog  node!  can  predict  system  response 
within  the  model  domain  of  three  degrees  of  freedom  without  knowledge  of  the 
numerous  and  varying  trunk  and  cushion  parameters  if  th^  mode  spring  and  damper 
coefficients  are  provided.  the  model  response  correlates  with  Jir.divik  test  data 
much  tetter  in  pitch  and  roll  than  in  leave.  Doth  the  pitch  and  roll  node  spring 
and  damper  coefficients  are  shown  to  be  linear  in  accordance  with  nod- 1  assumptions 
Model  response  in  heave  is  heavily  dependent  upon  the  mode  damping  ratio  (  X  )  arid 
the  undamped  natural  frequency  (  Ur,  ),  but  the  results  are  reasonable,  that  is, 
response  frequency  correlates  to  test  data  and  response  displacement  is  in  the 
proper  dn^v_«.ion  tut  sometimes  ol  i.ie.-rtect  magnitude,  lne  heave  node  spring 
coefficient  is  shown  to  hold  to  the  linear  approximation.  I  he  inability  to  match 
magnitudes  exactly  suggests  that  tie  heave  damper  coefficient  is  nun-linear. 
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Preface 


This  report  Is  about  a  Vertical  Energy  Absorption  Model  (VEAM)  for 
an  Air  Cushion  Landing  System  (ACLS)  chat  was  developed  from  full-scale 
experimental  tests  of  an  Australian  Jindlvik  drone  aircraft*  The  tests 
were  conducted  by  the  Air  Force  Flight  Dynamics  Laboratory,  Mechanical 
Branch,  and  selected  test  cases  were  used  to  verify  the  predictive 
capability  of  the  model  in  its  three  degrees  of  freedom*  heave,  pitch 
and  roll.  It  is  my  hope  that  the  results  of  this  effort  will  find 
practical  use  in  the  Flight  Dynamics  Laboratory  and  will  help  in  the 
advancement  of  Air  Cushion  Technology. 

I  wish  to  express  my  gratitude  to  ray  academic  advisor.  Captain 
James  T.  Karam,  Jr.  for  the  finesse  with  which  he  managed  the  "care  and 
feeding"  of  the  student  mind  during  the  experience.  His  assistance  and 
guidance  were  invaluable. 

I  am  particularly  grateful  to  Major  John  C.  Vaughn  III,  Principal 
Scientist  of  the  Mechanical  Branch.  AFFDL--  for  his  assistance  throughout 
the  period  of  this  study.  A  special  note  of  appreciation  is  due  to 
Mr.  James  T.  Steiger,  also  of  the  Mechanical  Branch,  who  provided  much 
motivational  and  technical  assistance.  Without  their  aid  this  study 
could  not  have  been  accomplished. 

My  children,  Scott,  Kimberly,  and  Tiffany,  will  never  know  what 
they  contributed  toward  the  completion  of  this  study,  and  for  that  I 
wish  to  thank  them.  My  wife,  Dorothy,  knows,  and  with  that  knowledge 
she  remained  a  fortress  of  cadence  and  a  fountain  of  understanding. 

She  is  a  Gpecial  person,  and  I  want  very  much  to  thank  her... for  typing 
this  thesis. ..for  everything. 
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A> (tract 

A  mechanical  analog  Vertical  Energy  Absorption  Model  (VEAM)  le 
developed  to  predict  the  dynamics  of  an  Air  Cushion  Landing  System  (ACLS) 
in  the  vertical  dimension.  Three  degrees  of  freedom  and  thus  three 
primary  modes  of  oscillation  are  investigated!  heave,  pitch  and  roll. 

Data  from  Air  Force  Flight  Dynamics  Laboratory  Tests  of  a  full-scale 
Australian  Jindivik  drone  are  used  to  develop  and  verify  the  model.  As 
part  of  the  comparative  analysis  of  Jindivik  data,  a  computer  program 
is  developed  that  will  analyze  the  peak  data  of  an  underdamped,  second 
order,  sinusoidal,  response  and  evaluate  the  mode  response  characteristics, 
damping  ratio  ( *f  ),  undamped  natural  frequency  ( cJn.  )  as  well  as  the  mode 
spring  and  damper  coefficients.  Another  computer  program  is  developed  to 
analyze  the  same  response  characteristics  using  the  method  of  peak 
overshoot  and  peak  times. 

The  VEAM  study  demonstrates  that  the  use  of  a  mechanical  analog 
prediction  scheme  for  an  AIl—  Cushion  Landing  System  has  sufficient 
merit  to  warrant  further  investigation.  It  shows  that  a  mechanical 
analog  model  can  predict  system  response  within  the  model  aomain  of  three 
degrees  of  freedom  without  knowledge  of  the  numerous  and  varying  trunk 
and  cushion  parameters  if  the  mode  spring  and  damper  coefficients  are 
provided.  The  model  response  correlates  with  Jindivik  test  data  much 
better  in  pitch  and  roll  than  in  heave.  Both  the  pitch  and  roll  mode 
•prlng  and  damper  coefficients  are  shown  to  be  linear  in  accordance  with 
model  assumptions.  Model  response  in  heave  is  heavily  dependent  upon 
the  mode  damping  ratio  ( 'f  )  and  the  undamped  natural  frequency  (tdn  ), 
but  the  results  arc  reasonable j  that  i9,  response  frequency  correlates 
to  test  data  and  response  displacement  is  in  the  proper  direction  but 
•ometlmcs  of  incorrect  magnitude.  The  heave  mode  spring  coefficient  Is 
shown  to  hold  to  the  linear  approximation.  The  inability  to  natch 
magnitudes  exactly  suggests  that  the  heave  damper  coefficient  is  non-linear. 
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AIR  CUSHION  LANDING  SYSTEM, 
DROP  DYNAMICS  THEORY  (MECHANICAL) 


I.  Introduction 


Air  Cushion  Development 

The  Air  Cushion  Landing  System  (ACLS)  is  a  natural  extension  of  the 
.  already  successful  application  of  air  cushion  theory  to  Air  Cushion 
Vehicles  (ACV)  and  Surface  Effect  Ships  (SES).  Bell  Aerospace,  a 
Division  of  Textron,  Inc.  and  a  pioneer  in  this  field,  in  concert  with 
the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL),  has  developed  and 
flight  tested  an  ACLS  on  the  I.ake  LA-4  aircraft  (Figure  (1)  as  early  as 

August  1967-  These 
tests  proved  that 
an  aircraft  could 
takeoff,  land,  and 
ground  maneuver  on 
a  cushion  of  air 
much  in  the  same 
manner  as  the  ACV's. 
Further,  due  to  the 
large  surface 

"contact"  area  of  an  ACLS,  vehicle  weight  and  landing  loads  were  widely 


1  ^  ►  •*  4  K  A  <1  A,,l  ^  «  M  M  4  M  1  M  •  4?  M  Ate  w  «  ■>  te  M.  A  A  A  A  •  •  A.  MM  {  M  _  m  — .  m  .  »  1  A.  I.  te  A  ^  f\  P  ^ 

Vi  J.O  wa.  A.  Dutcu  p  i.11  iuw  i.  uuupt  i-u  V  Vd^plUALiiidLCi/  KJ »  J  *“■  J 

psi).  In  certain  cases  this  characteristic  of  the  ACLS  will  motivate 
abandoning  conventional  landing  gear  on  aircraft  for  the  air  cushion. 

The  lew  footprint  pressures  free  the  ACLS  equipped  aircraft  from  a  hard 
surface  rur.vay  requirement  and  allows  routine  operations  on  unprepared 
surfaces,  marsh  land,  water,  and  even  obstacle  strewn  surfaces.  In 
cooperation  with  the  Canadian  government,  Bell  Aerospace  is  currently 
testing  an  enlarged  version  of  the  LA-4  ACLS  on  a  dt  .villand  CC-115, 
Buffalo  (Figure  2).  Success  of  the  Buffalo  tests  is  expected  to  result 
in  an  ACLS  test  series  on  the  Lockheed  C-130  Hercules. 
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Fig,  2,  The  deHavilUnd  CC-115,  Buffalo,  with  ACLS 

Air  Cushion  Concent  and  Definitions 

The  operation  of  an  ACLS  is  not  the  same  as  the  ACV.  The  uniqueness 
of  the  ACLS  requires  development  of  its  own  discipline. 

The  design  and  method  of  construction  is  in  no  way  fixed.  Many 
configurations  are  under  investigation  (Ref  3i26),  but  the  one  shown  in 
Figure  3  seems  to  be  relatively  standard.  The  shape  and  material  of  the 
trunk  has  a  significant  effect  on  its  capability. 

The  trunk  is  that  peripheral  portion  of  the  system  that  contains 
the  air  cushion.  Usually  an  elongated  donut  shaped  manifold  of  a  nylon/ 
rubber  material,  it  contains  the  flow  of  air  around  and  into  tne  cushion. 

The  cushion  is  comprised  of  that  volume  of  air  inside  the  trunk 
periphery |  that  is,  the  hole  of  the  donut. 

r» k  — .n  _  k.  .  u  »  j ..jj  .  j - * -  j  e  * j  _  . *.« 
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cushion  and  a  back-pressure  relief  valve  in  the  trunk  to  reduce 
springiness  on  landing.  The  trunk  and  the  cushion  together  absorb  the 
vertical  energy  of  a  descending  mass. 
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Early  model  and  inexpensive  trunks  are  made  of  a  non-stretchable 
material  that  functions  like  an  inflatable  beach  ball  which  maintains 
the  same  surface  area  whether  or  not  inflated.  More  sophisticated 
trunks  are  made  of  a  stretchable  material.  A  trunk  cannot  be  allowed 
complete  freedom  of  stretch  or,  like  a  balloon,  it  will  expand  under 
pressure  until  fracture.  An  ACLS  trunk  must  have  limits,  like  a  woman*# 
girdle,  which  allows  elastic  expansion  until  it  reaches  the  nylon  thread 
limit.  Unlike  a  girdle,  which  has  a  50l  stretch  capability,  a  Bell 
Aerospace  trunk  will  stretch  up  to  300Z  for  its  design  shape  limit.  This 
kind  of  trunk  can  be  installed  under  tension,  flush  to  the  aircraft 
undersurface,  eliminating  the  need  for  complicated,  heavy  storage  doors 
vhile  maintaining  a  low  drag  profile  when  not  in  use, 

Air  Cushion  Application 

The  Air  Force,  through  the  AFFDL,  is  showing  great  interest  in 
applying  air  cushion  technology  to  drone  recovery.  Most  drones  are  in 
the  same  size  and  weight  class  as  the  LA-4,  so  feasibility  has  already 
been  demonstrated.  Present  drone  operations  include  surface  launching 
from  a  rail,  dolly,  or  skid,  or  air  launching  from  under  the  wing  of  a 
mothership  (C-130).  Recovery  methods  include  mid-air  retrieval  by 
helicopter  or  parachute  descent  to  the  surface.  Costs  of  support 
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operations  and  limitations  on  the  numbers  of  drones  that  can  be 
handled  simultaneously  demand  a  better  method.  One  serious  proposal 
is  the  Air  Cushion  Recovery  System  (ACRS)*  A  recovery  system  differs 
from  a  landing  system  in  the  sens*  that  it  is  designed  for  landings  only 
and  does  nor.  need  the  higher  power  (airflow)  required  for  taxi  and 
takeoff  operations.  It  can  be  constructed  of  less  expensive  n on- 
stretchable  materials  which  can  be  stored  in  parachute  fashion  ready 
for  deployment  at  recovery.  In  a  current  test  series,  the  AFFDL  Is 
actively  engaged  in  the  development  of  an  ACRS  on  an  Australian 
Jindivik  drone  (Figure  4). 


Fig.  4,  The  Australian  Jindivik  Drone  with  ACRS 


Background 

The  technology  of  the  ACLS  is  in  its  infancy.  Major  research 
efforts  were  the  LA-4  program,  the  one  quarter  scale  and  the  one  tenth 
models  of  the  CC-115,  and  the  current,  full  scale  CC-115  program. 
Throughout  these  efforts,  the  primary  approach  to  development  of  a 
Vertical  Energy  Absorption  Model  (VEAM)  has  been  with  the  emphasis  on 
mass  air  flow  theory  (Ref  9).  This  approach  was  severely  complicated 
by  the  number  of  variables  changing  simultaneously  and  the  lack  of 
concrete  information  on  many  of  the  variables.  It  was  the-  initial 
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Jindivik  tests  that  suggested  an  alternate  approach)  drop  test  data 
Indicated  that  the  Jindivik  with  ACRS  was  behaving  like  an  underdamped, 
second  order  system.  In  April  1973,  the  AFFDL  decided  in  favor  of 
investigating  a  mechanical  analog  VEAM. 

The  mechanical  analog  promises  a  fast  reacting  prediction  scheme 
which  is  needed  as  a  computer  subroutine  to  be  joined  to  a  larger 
aircraft  simulation  program  that  is  being  developed  by  the  AFFDL  Flight 
Controls  Division.  The  ACLS  portion  of  the  program,  which  must  eventually 
include  forward  motion  with  its  associated  aerodynamic  and  braking  forces. 
Is  to  predict  the  forces  that  the  ACLS  contributes  to  the  overall 
dynamics  of  the  vehicle  on  landing.  A  real-time  dynamic  analysis  is  most 
important  in  order  to  keep  pace  with  other  force  inputs,  such  as  pilot 
stick  inputs,  and  to  maintain  a  realistic  simulation  of  vehicle  response. 

Purpose 

The  purpose  of  this  study,  under  the  sponsorship  of  the  AFFDL, 
Mechanical  Branch,  is  to  develop  a  Vertical  Energy  Absorption  Model  using 
a  mechanical  system  analog  that  will  successfully  predict  the  dynamics 
of  an  Air  Cushion  Landing  System  in  the  vettical  dimension. 

The  investigation  is  to  show  that  a  mechanical  analog  will  predict 
ACLS  response  even  though  the  physical  parameters  of  the  trunk  and  cushion 
are  not  known.  If  a  mechanical  analog  is  successful,  and  if  the  spring 
and  damper  coefficients  of  the  mechanical  system  eventually  can  be  related 
to  the  trunk  and  cushion  parameters,  then  this  method  of  analysis  promises 
a  simplified  design  approach  for  the  ACLS  without  the  need  to  first  design 
and  build  expensive  test  vehicles. 


Scope 

The  mechanical  analog  was  developed  solely  from  the  full  scale 
Jindivik  tests,  and  the  model  results  are  limited  by  the  nature  of  those 
tests.  Results  are  confined  to  cases  that  are  restricted  to  three  degrees 
of  freed<xnj  specifically,  heave,  pitch  and  roll.  This  feature  is 
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VEAM.  Not  an  explicit  part  of  the  analysis  is  the  understanding  that 
the  model  was  developed  from  a  non-stretchable  trunk  of  unique  design* 
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Extension  of  the  model  to  cases  where  a  stretchable  trunk  Is  employed 
Must  be  done  with  caution  since  the  energy  absorption  capabilities  and 
reactions  of  the  two  materials  may  vary  significantly*  The  basic  design 
of  the  Jlndlvik  trunk  differs  from  the  only  other  trunks  tested  to  date 
(the  LA-4  and  CC-115),  in  that  the  cushion  voluae  is  smaller  (in 
proportion  to  trunk  voluoe)  than  other  trunks*  The  full  effect  of  these 
differences  will  only  bacons  clear  after  more  tests  with  different 
configurations* 
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II.  Modal  Development  and  Theory 


The  Jindivik  Drone 

Tests  on  a  full-scale  Air  Cushion  Recovery  System  (ACRS)  were 
begun  by  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  in  February 
1973  at  Wright-Patterson  Air  Force  Base.  An  Australian  drone  aircraft, 
the  Jindivik,  was  provided  for  the  tests  and  fitted  with  a  Sandaire 
designed  ACRS  modified  by  the  AFFDL.  These  tests  motivated  the 
selection  of  the  Jindivik  drone  as  the  basis  for  the  Vertical  Energy 
Absorption  Model  (VEAM). 

The  Jindivik  was  set  in  a  test  cell  as  shown  in  Figure  S  and 
modified  slightly  to  accomodate  the  tests.  The  forward  hood  was 
removed  to  facilitate  instrumentation  and  the  landing  skid  was 
removed  for  the  ACRS  installation.  Weights  were  added  to  simulate  the 
Internal  equipment,  engine,  and  fusl  with  care  exercised  so  that  the 
total  weight  and  moments  of  inertia  of  the  test  vehicle  matched  its 
operational  counterpart.  Friraary  control  of  the  Jindivik  for  the  tests 
was  achieved  through  the  single  suspension  point,  chain  and  pully 
assembly,  attached  directly  above  the  center  of  gravity. 

A  planform  drawing  of  the  drone  in  Figure  6  locates  the  recording 
Instrumentation  used  in  the  AFFDL  tests.  Signals  from  this  instrumenta¬ 
tion  eventually  were  recorded  on  light-sensitive  paper  by  a  Honeywell 
vieicorder  (Model  906  A)  from  which  data  was  reduced.  Pressures  of  the 
trunk  and  cushion  and  the  displacement  of  the  nose,  wingtip,  and  center 
of  gravity  were  measured  on  all  tests  used  for  development  of  the  VEAM 
nodal. 

Changes  to  the  ACRS  trunk  design  component  fabrication,  and 
installation  were  the  cooperative  effort  of  the  AFFDL,  Centro  and 
Goodrich.  The  trunk,  depicted  in  Figure  7,  was  constructed  from  a 
non-stretchable  neoprene  coated  nylon  of  various  weights.  The  lewer 
surface  of  the  trunk  consisted  entirely  of  brake  tread  made  from 
three-eighths  inch  tire  tread  material.  The  heavy  brake  tread  results  in 
s  weight  penalty  (one-half  the  weight  of  the  entire  trunk  assembly),  but 
it  did  facilitate  precise  angular  positioning  of  vh'  ttunk  uuz*lca  for 
optimum  efficiency.  Two  Tech  Development  TD-530  ejectors  were  provided 
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e  -  Center  of  gravity  (CC) 

£3  -  Linear  displacement  transducer 
Q  -  Accelerometers 
Q  -  Centroid  of  ACPS 


Fig.  6.  The  Jindivik  with  ACHS,  Dimensions,  and  Instrumentation 
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Fig.  7.  The  Jindivik  ACRS  with  Dimensions 

for  the  portion  of  the  tests  used  to  develop  the  VEAM*  cue  provided 
mass  flow  to  the  trunk  and  the  other  to  the  cushion. 


Test  Method  and  Data 

Wine  AFFDL  Jindivik  tests  veto  used  to  develop  the  VEAM.  The  testa 
vers  In  the  three  primary  cedes  cf  oscillation  of  the  node  1 }  heave ,  pitch 
and  roll.  In  pitch  and  roll,  the  analysis  centered  on  perturbation  about 
the  hover  equilibrium  condition,  whereas  in  the  heave  case,  the  analysis 
vas  of  Jindivik  oscillation  dynamics  after  being  released  from  a  height 
greater  than  hover  equilibrium.  The  method  of  collecting  data  in  the 
three  inodes  influenced  the  analysis  of  each  mode  and  the  resulting  model. 

Roll  Mode  -  One  wing-tip  was  lowered  to  the  surface  and  released  at 
time  zero  (to  )•  The  wing-tip  linear  transducers  recorded  vertical 
deflections  of  the  wing-tip  as  oscillations  damped  back  to  hover 
equilibrium.  A  maximum  roll  angle  was  of  the  oraer  of  6  degrees. 

Pitch  Mode  -  Resistance  to  pitch  perturbation  was  significant 
enough  that  the  Jindivik  could  not  be  displaced  easily  and  released* 
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*s  with  the  roll  case,  Oscillation*  were  Introduced  until  the  pitch 
perturbations  were  of  sufficient  magnitude.  The  linear  transducer* 
attached  to  the  nose  boom,  recorded  vertical  deflection*  of  the  nose 
boom  relative  to  hover  equilibrium,  A  maximum  pitch  angle  of 
approximately  7  degrees  was  experienced. 

Heave  Hode  -  The  magnitude  of  heave  stiffness  was  large  enough  to 
»ake  it  technically  impractical  to  conduct  perturbation  testing.  Thus* 
the  Jindivik  was  raised  to  a  height  and  released.  This  test  method 
required  separation  of  the  analysis  into  two  regions*  the  first,  the 
free-fall  from  release  until  "contact"  of  the  ACRS  with  the  surface-  and 
the  second*  the  dynamic  oscillations  of  the  ACRS  about  the  hover 
equilibrium  condition.  Since  the  ACRS  rides  on  a  cushion  of  air* 
"contact"  does  not  necessarily  imply  the  physical  meeting  of  the  trunk 
and  the  surface.  Thus,  "surface  contact"  has  been  defined  as  In  Ground 
Effec-  (IGE)j  that  is,  the  position  at  which  the  first  distinguishible 
rise  in  trunk  or  cushion  pressure  can  be  detected.  Experimental 
testing  shows  tnac  trunk  ana  cushion  pressure  rise  within  a  fraction  of  • 
second  of  one  another,  and  can  be  considered  simultaneous  for  the  purpose 
of  defining  IGE-  Heave  data  was  complicated  by  the  fact  that  the  heave 
linear  transducer  could  not  be  placed  at  the  center  of  gravity  due  to  the 
fuselage*  nor  could  it  be  placed  directly  above  the  center  of  gravity  due 
to  the  suspension  hoist.  The  sensor  offset  resulted  in  heave  data 
being  affected  by  pitch  and  required  analytical  adjustment. 

A  summary  of  the  AFFDL  tests  in  its  three  significant  modes  of 
oscillation  and  under  various  conditions  is  presented  in  Table  1. 


TABLE  I 

AFFDL  JINDIVIK  TF5T  NUMBERS 
used  for  the  VEAM  development 


Test 

Mode 

Vent  Opened 

?c  * 
cc  <0 

Vent  Closed 

Pc  *  unk 
rac  ■  0 

Flow  Added 
PCp=12.5  psi 
mc  ■ 

Flow  Added 
fcp° 30.4  psi 
mc  ■  nt 

ROLL 

56 

55 

58 

PITCH 

5? 

60 

61 

HFAVE* 

123 

122 

121 

*  Jindivik  balanced  so  that;  CG  is  approximately  over  the  trunk  CP. 
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Tests  were  conducted  in  still  ambient  air  in  an  enclosed  test  cell* 
and,  as  described  above,  were  limited  to  perturbations  or  vertical 
descent  (no  forward  velocity).  Thu*,  no  aerodynamic  lift  or  drag 
forces  nor  brake  drag  forcer,  were  present  to  complicate  force  or  moment 
analys is . 


The  Vertical  Energy  Absorption  Model  (VvAtt) 

The  model  choser.  for  analysis  is  depicted  in  Figure  8. 


Fig.  8.  Schematic  of  Vertical  Energy  Absorption  Model 


Normally,  consideration  of  an  object  in  three  dimensional  rpace  results 
in  six  degrees  of  freedom  (DOF).  Jindivik  testing  methods  allowed  a 
reduction  to  T  JjOF.  The  elimination  of  air  flow  over  the  wings  and 
forward  movement  of  the  model  during  the  initial  test  series  had  the 
effect  of  eliminating  all  forces  in  the  horizontal  (>;,y)  piane.  Conse¬ 
quently,  it  is  assumed  that  there  is  no  lateral  displacement,  and  that 
all  movement  is  confined  to  the  vertical  (Z)  axis.  Because  there  are  no 
lift  and  drag  forces  to  create  side  loads,  yav»  is  not  present.  Thus,  the 
dynamic/perturbation  study  is  limited  to  the  vertical  mode  (heave)  and 
the  pitch  ar*4  roll  modes.  The  axis  system  chosen,  similar  to  n 
conventional  Euler  system  with  Z  positive  down,  has  its  or'.?ii<  at  the 
center  of  gravity  and  is  fixed  in  space  with  the  x,y  plane  parallel  to 
the  surface  vhen  the  center  of  gravity  ic  at  the  hover  equiiittiun 
condition.  This  seeming  complication  to  the  axis  is  required  in  order 
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to  generalize  the  analysis  and  yet  allow  a  changing  equilibrium  height 
(relative  to  the  sur.'ace)  as  the  mass  flow  into  the  cushion  changes. 

Model  Assumptions 

Objervation  of  the  Jindivik  oscillations  in  its  three  modes  results 
in  data  which  gives  insight  into  the  character  of  the  ACRS.  One  of  the 
first  Jindivik  dynamic  tests  conducted  was  the  roll  test  shown  in  Figure  9. 
The  plot  of  peak  displacements  suggests  an  underdamped  6econd  order 
response,  which  motivated  the  selection  of  a  mechanical  analog  model  for 
the  ACRS.  The  period  (  T  )  is  approximately  constant  throughout,  which 
Bakes  the  damped  natural  frequency  (  au  )  constant.  A  semi-log  plot  of 
the  peak  amplitude  excursions  from  equilibrium  demonstrates  the  linearity 
of  the  exponential  decay.  A  complete  discussion  of  this  analysis  is 
c^n^ained  i  Chapter  IV, 

There  is  -  jingle  effective  spring  and  a  single  effective  damper  that 
results  from  ea;-r.  mode  analysis.  With  the  assumption  that  the  springs 
and  dampers  are  linear  and  time  invariant,  it  is  possible  to  construct 
spring/dampe*-  units  (s/d  units),  each  a  linear  spring  and  damper  in 
parallel,  and  distj-ioute  them  as  depicted  in  Figure  8,  The  rotational 
roll  spring  and  damper  and  the  rotational  pitch  spring  and  damper  are  each 
divided  into  translational  s/d  units  displaced  at  ar.  assumed  moment  arm. 
The  linear  heave  3pring  and  damper  are  represented  by  the  combined  effect 
of  the  four  peripheral  s/d  units  and  the  center  cushion  s/d  unit.  The 
constants  of  proportionality  for  the  Jindivik  are  determined  from  the 
experimental  data  (Appendix  A). 

As  is  usual  with  lumped  linear  models,  analysis  requires  the 
concentration  of  mass  at  the  center  of  gravity  (point-mass  assumption). 

The  plane  surf  ~e  which  distributes  the  s/d  units  at  the  proper  moment  arm 
is  assumed  massless.  The  center  of  gravity  can  be  varied  from  the  center 
of  pressure  (also  the  geometric  center:  of  the  trunk)  by  the  distance  6  . 

In  the  model,  the  center  of  gravity  and  the  center  of  pressure  are  in  the 
same  plane. 

Approach  to  Analysis 

Whan  the  i*CltS  is  Cr.it  of  Ground  Effect  (OCc),  it  is  inoperatlvej  that 
ia,  pressures  of  the  trunk  and  cushion  are  constant  and  no  forces  or 


Fig.  9.  Linearity  of  Exponential  Decay  in  Roll 
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moments  are  present.  As  the  ACRS  descends  in  free  fall  toward  the 
surface*  it  reaches  the  region  (IGE)  where  reaction  with  the  surface 
causes  a  pressure  rise  in  both  the  trunk  and  cushion. 

The  instant  the  IGE  condition  is  encountered  marks  time  zero  (to  ) 
for  activation  of  the  analytical  model.  This  displacement  corresponds  to 
the  spring-neutral  position  (Jf^  )  as  depicted  in  Figure  10. 


Fig.  10.  Schematic  of  Vertical  Displacements 


Any  further  downward  displacement  causes  reaction  forces.  Descent 
of  the  model  continues,  compressing  each  spring  (AZ  )  until  the  combined 
upward  force  of  the  spring/daraper  units  (s/d  units)  supports  the  weight 
of  the  model.  This  displacement  is  the  equilibrium  position  (  )  and 

corresponds  to  equilibrium  hover  for  the  Jindivik.  When  the  center  of 
gravity  is  displaced  from  the  center  of  pressure,  tl  e  Jindivik  longitudinal 
centerline  will  not  be  level  with  the  surface.  A  non-zero  equilibrium 
pitch  angle  (  ©e^.  )  will  result  so  that  for  each  s/d  unit 

AZ  e  (Moment  Arm)  sin  ©e<^  (1) 

The  equilibrium  pitch  angle  is  non-zero  due  to  the  combined  physical 
constraints  of  (1)  the  CG  of  the  vehicle  being  aft  of  the  trunk  CP  in 
normal  configuration.  (2)  the  assuanrinn  that  the  roll  and  the  pitch 
spring  coefficients  are  symmetrically  placet!  relative  to  the  trunk  CP,  and 
(3)  the  assumption  that  the  symmetrical  spring  coefficients  are  equal 
in  magnitude. 
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Linearity  of  Springs 

The  curves  in  Figure  11a  illustrate  that  a  linear  approximation  for 
a  non-linear  spring  is  usually  valid  within  a  small  perturbation  of  a 
given  position,  in  this  case,  equilibrium.  Figure  lib  illustrates  that 
real  spring  phenomena  is  not  infinitely  linear,  but  that  non-linearities 
are  visually  present  at  the  extremes.  For  the  VEAM,  when  it  is  above  the 
contact  level,  no  forces  should  exist.  At  this  extreme  the  spring 
relationship  F=  implies  that-&-»o  at  displacements  less  than  contact 
(Xc,  )•  At  the  other  extreme,  where  the  vehicle  underbelly  strikes  the 
surface,  the  spring  relationship  again  terminates.  It  is  suspected  that 
the  spring  goes  non-linear  prior  to  fuselage  and  surface  contact  (Xs), 
These  two  extremes  define  the  limits  of  the  domain  for  successful  VEAM 
operations. 


Fig.  11.  Spring  Linearity 


It 


III*  The  Model  Equations  of  Motion 


The  Force  Equation 

In  this  chapter,  the  equations  of  notion  that  apply  to  the  three 
degrees  of  freedom  will  be  developed.  Each  spring/damper  unit  M 

illustrated  in  Figure  12  will 
result  in  a  force  with  constants 
of  proportionality  for  displace¬ 
ment  Jit  and  velocity  -Ci  for  each 
mode.  The  total  force  described  as 

Fj  -  -  (  A  i  Zc  j  +  -<^i  Z-ci)  (2) 
where  displacement  is  measured 
positive  down  from  the  spring- 
neutral  or  contact  position* 

When  the  model  is  released  from 
the  spring  neutral  position,  which 
is  defined  as  horizontal  to  the 
surface,  it  will  settle  downward 
due  to  model  mass  and  seek  an 
equilibrium  position  so  that  the  forces  and  moments  for  each  degree  of 
freedom  ai i  exactly  balanced.  This  equilibrium  position  can  vary  in 
relation  to  the  horizontal  surface  depending  upon  the  relative  location 
of  the  s/d  units  to  the  center  of  gravity  and  the  center  of  gravity  to  the 
center  of  pressure  of  the  trunk.  Thus,  at  each  s/d  unit  attachment 
point  2  ,  if  the  variable  deflection  were  defined  as  AZ  ,  then  the  total 
force  at  each  point  would  be 
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Fig.  12.  A  Spring/Damper  Unit 
in  Arbitrary  Deflection 
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However,  for  any  given  model  configuration,  the  forces  and  moments  due  to 
the  deflection  A7.  from  spring-neutral  to  equilibrium  are  constant  so  as 
to  exactly  counteract  '.he  presence  of  the  mass  at  all  times.  If  the 
deflection  from  contact  lo  equilibrium  were  eliminated  from  an  analysis, 
the  resulting  equations  represent  a  dynamic  force  balance  relative  to 
the  equilibrium  condition  where  each  s/d  unit  force  contribution  would  be 
represented  by 


Fi 


-V  -Ci  )  , 


(4) 
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In  this  last  eauatfon,  represents  the  vertical  deflection  of  each 
s/d  unit  attachment  point  from  the  equilibrium  position  (  ”  Ze}{  )• 

Perturbation  from  Equilibrium 

Analysis  of  the  model  in  its  three  degrees  of  freedom  starts  by 

considering  the  forces  that  act  upon  it  in  time  (equation  A)*  The 

forces  remaining  after  model  weight  is  eliminated  are  dependent  upon  ^ 
« 

and  ^  •  The  model  is  given  a  positive  deflection  in  heave,  pitch  and 
roll  simultaneously.  Allowing  that  all  angles  are  small  (Appendix  C), 


—7'  ^  /— 


Fig,  13.  The  VE AM  in  Arbitrary  Deflection  (Heave, Pit'.h  and  Roll) 
Figure  13  illustrates  that  the  deflection  at  each  s/d  jntt  location  is 
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The  rate  of  deflection  of  each  s/d  unit  point  is  the  derivative  of 
equations  (5). 
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J,  =  z  -  (4  *  4)  e . 

L=i“^i4>“6e 

J>*z+(&-5)e 

c«) 

}«  =  ±  -  <S  6 


Substitution  of  equations  (5)  and  (6)  into  equation  (A)  provide  the 
necessary  information  for  a  force  and  moment  analysis. 


The  Heave  Mode 

Summing  all  forces  in  the  vertical  direction,  where  the  forces  are 
as  previously  described,  results  in 


IFz  = 


(7) 


Substitution  of  the  complete  s/d  unit  force  equations  (the  combination 
of  equations  (A),  (5)  and  (6),  into  equation  (7)  and  separating  like 
orders  of  heave  displacement  and  pitch  angle  results  in  the  heave  mode 
equation 


Z 


z  + 


-free  __ . 
m 


e 

rr\ 


ie 


6  e 


(8) 


where  the  effective  spring  and  damper  coefficients  are 


-Ce  - 

(9) 

A 

Mg  - 

The  Roll  Mode 

In  a  like  manner,  the  roll  equation  develops  from  a  summation  of 
moments  about  the  longitudinal  x-axis. 

<)>  =  -  F*(H)  .  oo) 

In  this  case,  the  physical  symmetry  of  the  trunk  about  the  longitudinal 

axis  suggests  that  the  corresponding  s/d  units  are  also  symmetric}  or 
•Ci  -  -O  h  =  Cj(h 

-fta.  —  Jl*  — 

The  complete  s/d  unit  force  equations  are  substituted  into  (10)  and 


19 


the  coefficients  of  like  order  derivatives  are  collected  to  jive 


•  $  +  <H 


-c».m  b*  -i 

7.  Tka  J 
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'o  1 
£  Tw-1 


O 


(12) 


The  Pitch  Mode 

Development  of  the  pitch  equation  Involves  a  bit  more  algebra,  but 
the  method  is  the  same.  Moments  are  summed  about  the  lateral  y-axis. 

*fsS  as) 

and  the  complete  s/d  unit  force  equations  are  substituted.  To  the 
assumption  of  the  roll  s/d  units  (equation  12)  is  added  the  more  tenuous 
assumption  that  the  fore-aft  s/d  units  are  also  symmetric,  or 


•Ci  -  -  *C- 1.1 

ill  -  Jfe*  *  ill,} 


(14) 


Equations  (14)  are  considered  more  assumptive  due  to  the  slightly 
different  volume  between  the  fore  and  aft  trunk  sections  (see  the  cross- 
sectional  areas  in  Figure  7),  but  the  difference  is  considered  negligible. 
The  complete  pitch  equation  is 


e  -3  e  ^  [z-e  J  (If*  £2]  -  *  -t^}  ] 

+  e>{4rl  zJUf&M'l  +2  *,.£‘4  iu/h 


Jyy 


iz 

Tyy 


Ce 


(15) 


Free  Fall  Analysis 

The  free  fall  case  requires  a  separate  analysis  in  a  different 
time  frame.  After  implementing  the  mode  assumptions  presented  in 
Chapter  II,  a  free  body  diagram  of  the  Jindivik  Out  of  Ground  Effect 
(OGE)  shows  that  weight  is  the  only  force  acting  on  the  system  as  it 
falls  in  the  vertical  dimension.  Time  is  measured  from  release  and, 
so  as  not  to  confuse  with  the  time  domain  for  the  model  ICE,  free  fall 
time  is  designated  T  .  The  height  of  the  vehicle  in  free  fall  is 
designated  ( Jh.  )  and  follows  the  terminology  presented  in  Figure  10. 
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It  is  measured  positive  down  to  maintain  the  same  sign  convention  as 
th«  model.  * 

The  equations  which  describe  free  fall  with  no  initial  velocity  are 


7  •  x  3  't'1  *♦*  Ai 

(a) 

1  -  3T- 

(b) 

2  «  -3 

(c) 

This  free  fall  analysis  is  used  in  conjunction  with  the  heave  tests 
as  previously  described,.  The  values  of  the  release  height  (  iui )  and  the 
time  to  initial  contact  (  T~c  )  are  determined  experimentally,  and  from 
(16a)  the  height  at  contact  (  JLC )  is  determined.  Equation  (16b)  provides 
another  initial  condition,  the  velocity  at  contact.  These  values,  with 
coordinate  transformation,  are  then  used  in  an  analysis  of  the  character¬ 
istics  of  the  equations  of  motion. 

Summary 

The  complete  set  of  equations  to  describe  the  model  in  its  three 
degrees  of  freedom  is  the  combination  of  equation  (8),  (12)  and  (15). 

It  should  be  noted  that  both  the  heave  equation  (8)  and  the  pitch 
equation  (15)  are  coupled  to  one  another  if  the  center  of  gravity  of  the 
vehicle  is  not  at  the  center  of  pressure  of  the  trunk,  ie,  £tO,  This 
ia  to  say  that  a  pure  heave  input  motion  will,  in  time,  result  in  a 
pitch  motion  that  may  remain  even  after  the  heave  notion  has  damped  out, 
and  vice-versa. 

Solution  of  these  equations  of  motion  is  the  subject  of  the  next 
chapter,  but  it  must  be  remembered  that  the  dynamic  analysis  is  in 
reference  to  the  equilibrium  position  of  the  vehicle.  The  equilibrium 
position  can  vary  according  to  the  vehicle  configuration  and/or  flow 
rates  into  the  trunk  and  cushion.  In  heave,  the  difference  in  height 
between  equilibrium  and  the  surface  is  constant  for  any  given  test. 
Program  DATANYL  in  Appendix  A  computes  this  value.  For  roll,  the 
equilibrium  axis  is  parallel  to  the  horizontal  axis  at  all  times  due 
to  the  fact  that  the  trunk  is  symmetrical.  Pitch  is  the  same  as  roll 
only  if  the  center  of  gravity  of  the  vehicle  is  coincident  with  the 
center  of  pressure  of  the  trunk.  For  most  pitch  cases,  where  the  CG  is 
behind  the  CP,  an  equilibrium  pitch  angle  (0^),  relative  to  the 
horizontal,  will  result.  An  expression  for  0e<^  is  in  Appendix  B. 
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IV.  Analysis  of  Jindivik  Data  (Theory) 


Analysis  of  Roll  and  Fitch  Data 


The  equations  of  notion  developed  in  Chapter  III  can  be  related  to 
a  standard  form  for  linear,  constant  coefficient,  second  order 
differential  equations  (Ref  It 244)  where  the  coefficient  of  the 
displacement  term  is  equal  to  and  the  coefficient  of  the  first 

derivative  of  displacement  is  equal  to  For  example,  in  the 

case  of  roll,  these  are  written  as 

-C  0**  (17) 

-9?  »,q  _  , 

The  parameters  u)*\  •  the  undamped  natural  frequency,  and  f  ,  the 
damping  ratio,  serve  to  characterize  the  differential  equation  and 
thus  the  system.  The  solutions  to  the  VEAh  differential  equations  of 
motion  (equations  8,12  and  15)  are  of  the  form 

PCt)  =  Cb  e~1u)r,t  sin  (wdl- +<$)  0«> 

where  DCt)  is  angular  or  linear  displacement  as  a  function  of  time,  Cb 
is  a  constant,  and  uW  is  the  damped  natural  frequency.  By  definition, 

0)4  =  1  -  -f*-'  <l9> 

The  nature  of  the  system  being  analyzed  depends,  to  a  large  part,  on  the 


Fig.  14.  Ceneral  L’nderdamped  Sinusoidal  Response 
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value  of  "f  .  For  the  case  of  this  study  where  f  <  ^.0  ,  the  system  la 
called  "underdamped"  and  the  natural  response  appears  as  in  Figure  14. 
The  period  of  system  response  (  T  )  is  the  time  required  to  traverse 
one  full  cycle,  thus  T  »tf1  -  tp<  -  1  / $&  ftps  .  From  this,  the 
following  relationship  is  established 

UM  (rad/sec)  *  2Vt  *  **  ^  (cyc/sec)  (20) 

The  period,  which  can  be  found  from  experimental  data,  leads  directly 
to  a  value  for  the  damped  natural  frequency. 

A  description  of  the  system  as  it  damps  to  the  equilibrium 
position  is  available  by  comparing  two  adjacent  peak  amplitudes*  This 
comparison  is  called  the  amplitude  ratio  and  defined  as  Di , 

Combining  equation  (19)  with  (18)  for  two  adjacent  peaks  and  solving 
for  the  amplitude  ratio  gives 

Oi  fj/lTFT*-' 

—  ^  e  <a> 

The  logarithmic  decrement  (  <i  )  is  equal  to  Lhe  natural  logarithm  oi  the 


Oh 


2  IT-  f 


amplitude  ratio,  and 

^  ^  ~oI  W=W 

Alternately,  the  damping  ratio  is 

i* 


(22) 


T  =  l 


f 


(lit)*  -  dZ 

The  preceeding  analysis  was  computerized  to  evaluate  the  system  in 
roll  and  pitch,  to  provide  the  mode  spring  and  damper  coefficients. 


(23.' 


and  to  plot  the  data. 
Appendix  A. 


This  program  and  its  results  are  presented  in 


Analysis  of  Heave  Data 

The  initial  analysis  of  Jindivik  heave  data  indicated  rapid 
damping  of  the  heaye  mode  with  another  mode  remaining.  This  is 
illustrated  in  Figure  15  .  Evaluat  ion  of  the  data  past  the  mode  change 
point  indicated  that  the  pitch  mode  characteristics  now  predominated. 
Test  results  ‘.tu  that  the  heave  mode  damped  out  in  about  three 

peaks  (P  .75  to  1.45  seconds).  With,  so  few  peak  amplitude  points,  the 
heave  data  was  considered  too  sparse  to  analyze  in  the  came  manner  as 
pitch  and  roll.  For  thi3  reason,  an  analysis  based  on  percent  overshoot 
and  rise  time  was  selected  to  characterize  the  system.  This  method  of 
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solution,  which  is  illustrated  in  Figure  16,  accounted  for  the  fact  that 
the  model  has  both  initial  displacement  and  initial  velocity  at  model 
activation  (surface  contact). 

Starting  with  the  heave  differential  equation  of  motion  (8), 

Laplace  transforms  are  used  to  find  the  solution 


Z<t)  -fwnt 

x —  -  C  (cos  u)dt  ♦  K  stntudt) 


(24) 


where 


f 

K  =  VT^f 1 


2e  L*n  ll-  t’*'  (25) 

By  using  a  trigonometric  double-angle  identity,  (24)  can  be  simplified 
to  the  form 

-  IvTk* 

CrO 

-1 


-  Tu)nt.  .  ./\ 

G  S\n  (um  t  +  q>/ 


where 


4>' 


CO  u 


H 


s»n 


(26) 

(27) 


Maximum  overshoot/undershoet  occurs  when  S\n  (u)d  fc  -  -  1-0. 


Thus, 


Z(t)rnax 


=  ill 


-f  CJn  t, 


7.  •*  *'  w  (28) 

Equation  (26)  is  differentiated  with  respect  to  time  and  solved  in 
order  to  obtain  a  general  expression  for  peak  time  (  tp )  v3iid  for  any 
response  peak  depending  on  the  value  of  the  integer  (r\), 

"t»p  —  Zq  /  i*)n  (Zof  "*■  ZokJn)  viTI"  /  1  ~  T2-'  (29) 

Equation  (28)  is  solved  for  the  damping  ratio  and  equation  (29)  for 
the  undamped  natural  frequency  to  provide  two  equations  suitable  for  an 
iterative  solution  on  a  computer. 

(30) 


,  -J?n.[|Z(tw*  /Zol  /VW-Kr] 

T  =  - 


tp 


(a) 


a  7.0  (  T  Zc  +  ZoL*>n) 


tJn 


tp  (fZe  +  ZciOrx') 

The  program  which  accomplishes  this,  called  RISE,  is  presented  in 
Appendix  A. 


(b) 


Results 

The  results  of  the  analysis  of  Jindivik  data  from  both  programs 
DATANYl  and  RISE  are  presented  in  the  following  tables.  Table  II 
characterizes  the  system  for  each  node  for  each  AFFDI.  rest,  and 
Table  III  summerizes  results  that  are  peculiar  to  the  heave  mode  only. 
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TABLE  II 


KESHA*  OF  ANALYSIS  OF  JINOIVIK  DATA 


Program  DATANYL 


Mode  Case  Test 
Ho 


Pitch  I  1  59 


Roll 


Program  RISE 


123AV 
122AV 
1 21 AV 


Zeta 


Wn 

Mode  Coefficient 

(rad/sec) 

(ft-lbf/rad) 

cE 

!  (ft-lbf/rad/sec) 

6.66 

80,394.76 

1,052.36 

6.20 

69,524.92 

1,036.00 

5.66 

57 , 917 . 22 

731.38 

1.84 

4,040.29 

234.06 

1.64 

3,077  99 

221.05 

1.52 

2,753.88 

227.53 

Heave 

1 

123 

3 

122 

4 

121 

First  Peak 


.3482  15.22 
.3533  16.61 
.3771  15.48 


Average 


14.90 

14.67 

14.24 


Cases 

1  -  Vent  opened 

2  **  Vent  closed 

3  -  Flow  added  (mc 

4  -  Flow  added  (mc 


Clbf/ft 


17  >769.30 
21,163.15 
1C, 381. 59 


17,029.96 

16,506.26 

15,554.68 


(lbf-sec/ft) 


813.04 

900.29 

895.57 


597.99 

410.06 

393.02 
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V.  Anal/sis  of  the  Spring/Damper  Units  (Theory) 


General 

The  analysis  of  the  model  equations  of  notion  is  completed  by 
computer  solution.  The  computer  analysis  is  presented  in  Appendix  B. 
Prior  to  implementing  this  solution,  the  spring/damper  units  must  be 
fully  described. 


SprinR/Pamper  Unit  Evaluation 

If  it  were  not  for  the  coupling  terms  in  the  heave  and  pitch  cases, 
the  equations  of  motion  could  be  solved  with  the  mode  characteristics 
only  (f  *s  and  uJn*s).  The  coupling  requires  knowledge  cf  the  individual 
s/d  units. 

In  Chapter  11,  it  was  emphasised  that  only  one  "spring’Value  and 
one  "damper"  value  could  be  derived  from  the  data  of  each  node.  Analysis 
of  two  of  the  modes,  pitch  and  roll,  resulted  in  torsional  springs  and 
dampers,  each  of  which  now  can  be  transformed  into  a  pair  of  linear 
springs  and  dampers  at  an  assumed  moment  arm  to  satisfy  the  model 
configuration.  The  moment  arm  chosen  for  this  analysis  is  based  on  the 
Jindivik  ACRS  trunk  (Figure  7).  The  tangential  contact  points  of  the 
trunk  and  the  surface  with  the  xz  and  yz  planes  define  the  points  for 
measurement  of  trunk  length  and  width.  Assuming  trunk  symmetry,  the 
proper  moment  arms  are  half  of  the  length  (-wl.)  and  half  the  width  (b/2.). 

In  Chapter  IV.  the  relations  between  the  mode  characteristics  and 
the  cpring/danper  coefficients  were  established.  The  example  for  the 
roll  case  is  repeated  here. 


Jfc  -  [  T«  uw2-  J  /  t*- 


00 

(17) 

(b) 


Expressions  for  the  heave  and  pitch  cases  can  be  found  by  similarly 
equating  the  coefficients  of  the  differential  equation,  but  the  heave 
and  pitch  spring/damper  coefficients  are  functions  of  each  other. 
Simultaneous  solution  of  the  expressions  result  in  a  form  that  is  in 
terms  of  the  mode  characteristics  only.  Thesa  arei 
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A?  - 
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(32) 

(33) 

(34) 


Numerical  solution  of  the  s/d  units  is  accomplished  in  program 
MODANYL  (Appendix  B),  the  results  of  which  are  presented  in  Table  IV. 


TABLE  IV 


LINEAR  SPRING/DAMPER  UNIT  COEFFICIENTS 
(from  Program  MODANYL) 


Mode-Test  No. 

Linear  s/d  Units 

■f  and  u)n 

Position 

Spring  Coefficients  (K) 

Damper  Coefficients 

(lbf/ft) 

(C) 

(lbf/ft/see) 

H-122 

1.3 

1,174.8 

6.9 

P-61 

2,4 

639.1 

52.9 

R-58 

5 

17,340.1 

772.4 

H-122AV 

1,3 

1,232.9 

13.0 

1  P-61 

2,4 

639.1 

52.9 

R-58 

5 

12,623.1 

274.6 

H-123 

1,3 

1,804.9 

16.7 

P-59 

2,4 

936,5 

54*4 

j  R-56 

5 

9,434.6 

599.5 

i 

j  H-123AV 

1,3 

1,780.2 

18.6 

P-59 

2,4 

936.5 

54.3 

R-56 

S 

11,439.4 

466.5 

Note  !«  For  the  heave  case,  calculations  are  based  on  data  acquired 
where  the  CG  was  shifted  above  the  trunk  CP. 

Note  2»  For  the  heave  rr.ode;  the  first  peak  analysis  was  used  for 

and  Da  unless  subscripted  by  AV,  in  which  case  the  average 
analysis  was  used. 

Note  3i  Linear  s/d  unit  values  are  based  on  the  Jindivik  Trunk, 

Figure  7,  where  £/2.  a-  53.85  inches  and  b//L  ®  17.6  inches. 
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VI.  Analysis  of  the  Vertical  Energy  Absorption  Model 


General 

The  Vertical  Energy  Absorption  Model  vis  exercised  on  a  CDC  6600 
computer  with  the  goal  to  duplicate  the  Jindivik  test  of  Table  I.  The 
model  results,  which  compare  favorably  with  the  original  tests  and 
with  one  another,  are  presented  in  graphic  form  in  Appendix  B  with  the 
exception  that  the  responses  from  a  single  test  condition  (flow  added- 
tests  58,  61,  and  123)  are  included  in  this  chapter  in  order  to  analyze 
them  more  closely.  A  variable  is  present  in  these  comparisons  since  the 
heave  tests  were  conducted  with  the  CG  approximately  above  the  center  of 
pressure  (CP)  of  the  trunk,  where  as  the  roll  and  pitch  tests  were  in 
normal  configuration  with  the  CG  aft  of  the  CP.  The  exact  effect  of  the 
change  of  the  CG  location  on  the  heave  data  is  not  known. 


Spring  Coefficient  Comparison 

Before  examining  the  model  responses,  it  is  worthwhile  to  compare 
the  mode  values  of  the  spring  coefficients  from  Table  II  with  earlier 
experimental  AFFDL  static  test  results  which  are  summarized  in  Appendix  D. 
These  values  are  plotted  in  Figure  17  in  a  kind  of  histogram  where  the 
abscissa  represents  the  point  cushion  conditions  of  vent  opened  (VO),  vent 
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Fig.  17.  Comparison  of  opring  Coefficients 
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closed  (VC)  and  flow  added  (FA).  Although  point  conditions*  the 
abscissa  does  suggest  a  scale  of  increasing  cushion  effectiveness  (  Pc  and 
>Tv.o  increasing).  This  suggests  that  the  oynamic  roll  spring  coefficient 
and  the  static  roll  stiffness  are  comparable*  Thus*  for  an  approximation* 
the  dynamic  roll  spring  coefficient  can  be  found  from  the  easier  and  less 
expensive  static  tests.  Another  advantage  is  that  the  static  analytic 
techniques  are  easier  than  the  dynamic  ones.  The  dynamic  pitch  stiffness 
also  shows  reasonable  agreement  between  the  dynamic  and  static  values  and 
does  approach  the  static  value. 

The  heave  case  shows  dissimilar  results.  The  first  peak  analysis 
values  appear  to  increase  while  the  average  analysis  values  show  the 
same  decreasing  trend  with  increasing  mass  flow  as  the  static  test  values. 
In  both  cases*  there  is  a  definite  magnitude  difference  in  the  values 
between  the  dynamic  and  static  values  of  the  spring  coefficients  for  any 
given  test  configuration.  This  fact  suggests  the  presence  of  a  dynamic 
spring  phenomenon  that  was  net  observable  during  the  static  testing: 

Model  Roll  Response 

The  test  values  of  the  damping  ratio  *f  and  undamped  natural 
frequency  u)n  used  to  characterize  model  response  are  those  in 
Table  II.  In  all  roll  tests,  *f  generally  increases  and  uVi  generally 
decreases  with  increasing  mass  flow  to  the  cushion. 

In  Figure  18  for  roll  test  58,  the  response  of  the  model  shows 
excellent  correlation  to  the  test  data  in  all  respects. 

Model  Pitch  Response 

For  the  pitch  tests.  Table  II,  "f  *s  show  no  decisive  trend,  while 
tOn’s,  like  the  roll  mode,  show  a  decreasing  trend  with  increasing  mass 
flow  to  the  cushion. 

Figure  19,  for  pitch  test  61,  generally  shows  good  results  except 
it  appears  that  the  equilibrium  axis  of  the  response  is  angled  slightly 
to  the  equilibrium  axis  of  the  data.  Referring  momentarily  to  Figure  30, 
page  60  ,  in  Appendix  A,  the  equilibrium  plots  <•>  of  the  original  data 
show  a  definite  increasing  trend  with  time.  The  equilibrium  height  (HEQ) 
of  each  test  was  defined  as  the  average  value  of  these  individual 
equilibrium  positions.  Back  to  Figure  19,  it  is  the  average  equilibrium 
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Fig.  IF..  Model  Response,  Roll  Teat  58,  Wing  Tip  Transducer 
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height f  a  constant,  that  was  used  in  program  MODANYL  to  adjust  the 
original  data  to  the  same  equilibrium  axis  as  the  model  response  so  that 
the  two  could  be  compared.  Thus,  a  slight  deviation  is  introduced  to  the 
plot.  The  magnitude  of  this  deviation  is  equal  to  the  difference  between 
the  average  equilibrium  height  and  the  actual  equilibrium  height  at  an7 
time* 

Model  Heave  Response 

The  heave  response  requires  several  computer  runs  with  different 
initial  conditions  in  order  to  evaluate.  In  Chapter  III,  the  axis 
chosen  for  the  development  of  the  equations  of  motion  was  based  on  the 
equilibrium  position.  Experimental  data  showed  that  the  equilibrium 
longitudinal  axis  made  an  angle  of  about  3-4.5  degrees  nose  up  from  the 
horizontal  plane  when  the  CG  was  displaced  aft  of  the  CP.  When  the 
aircraft  Is  dropped  level  with  the  surface,  then  an  additional  initial 
condition  in  pitch  is  established  where  ©  ~  -  3.0  degrees  relative  to 
the  equilibrium  axis  of  the  model.  The  first  peak  and  average  f's  and 
uij's  used  to  examine  the  heave  mode  response  are  those  calculated  in 
program  RISE  (Appendix  A)  and  presented  in  Table  V. 


TABLE  V 


7£TA  &  un  COMPUTATIONS  FOP.  HEAVE 
(Program  RISE,  modified) 


Peak 

TEST 

123 

TEST  122 

TEST  121 

ZETA 

Wn 

ZETA 

Wn 

ZETA 

Wn 

1 

.3482 

15.22 

.3533 

16.61 

.3771 

15.48 

2 

.2444 

14.72 

.1627 

14.40 

.1691 

14.52 

3 

.1922 

14.77 

.1296 

14,35 

.0851 

13.93 

4 

.1575 

13.80 

.1078 

13.32 

5 

.1847 

14.77 

.2515 

14.37 

6 

.1057 

14.12 

.0390 

13.82 

Average 

Valuer. 

.2616 

14.90 

.1822 

14.675 

.1799 

14.24 
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In  the  first  two  heave  responses  presented,  the  model  is  exercised 
with  initial  conditions  on  heave  only  (  4>  ■  a  ■  O)  and  thus  no  pitch 
coupling  develops.  The  results  of  these  initial  conditions  are  what  was 
expected  from  che  Jindlvik  tests  when  the  CG  was  moved  to  a  position 
above  the  trunk  CP.  The  fact  that  the  Jindlvik  test  data  shows  tha 
heave  test  coupling  into  a  pitch  mode  suggests  that,  during  the  actual 
tests,  the  CG  was  not  exactly  over  the  CP  and/or  the  initial  pitch 
condition  was  not  exactly  level  at  contact. 

In  Figure  20a,  the  ~f  and  from  the  first  peak  analysis  are 
used.  The  response  correlates  well  with  the  first  data  point  but  is 
too  heavily  damped  to  approach  the  other  data  points. 

In  Figure  20b,  the  average  values  for  -f  and  u;*  are  used.  Even 
though  the  response  overshoots  the  first  peak,  a  much  better  subsequent 
response  is  achieved. 

In  the  next  two  cases,  suspected  actual  Jindivik  test  conditions 
are  duplicated;  that  is,  the  model  is  activated  parallel  to  the  surface, 
which  puts  an  initial  condition  on  pitch,  ar.d  with  initial  conditions 
on  heave.  In  these  response  examples,  the  model  resembles  normal 
Jindivik  configuration  where  the  CG  is  behind  the  CP. 

In  Figure  20c,  the  first  peak  analysis  shows  a  degradation  of 
response  to  the  first  peak  due  to  pitch  coupling  and  a  tendency  for  the 
pitch  mode  to  dominate  the  response)  that  is,  heave  is  not  effective 
after  the  first  peak.  It  is  obvious  that  the  pitch  peaks  are  nowhere 
near  the  data  points,  but  this  tiiy  not  be  too  serious.  The  period  of  the 
pitch  mode  is  constant  and  approximates  the  period  of  the  data.  This 
means  that  the  damped  natural  frequences  ( Ld<l  )  do  compare.  In  addition, 
the  peak  amplitudes  of  the  pitch  response  fit  exactly  into  the  data 
exponential  decay  curve,  thus  suggesting  a  comparable  logarithmic 
decrement  and  pitch  damping  ratio  ). 

In  Figure  20d,  the  same  initial  conditions  as  Figure  20c  are  used 
but  with  the  average  values  for  "f  andu;.,.  The  depredation  of  the  first 
peak  response  actually  assists  toward  matching  the  data  point.  Subsequent 
response  is  much  improved  in  heave,  with  the  heave  node  matching  or, 
at  least,  moving  toward  all  data  points. 

A  closer  inspection  of  the  character  of  the  heave  mode  is  possible 
with  Table  V  which  incluies  a  summary  of  -f  and  calculations  from 
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Fig.  20a.  Modal  Response,  Ueava  Test  122,  CC  Transducer 
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program  RISE  for  all  peaks*  In  all  the  Leave  test  cases,  the  first  peak 
values  of  -f  an',  u)n  are  significantly  I’tgher  tlian  the  subsequent  values. 
The  subsequent  values  of  cJn.  are  essentially  ccnscamt.  Equations  (l?b), 
(32)  and  (34)  imply  that  when  iAi  is  constant,  then  the  spring  coefficient 
is  constant,  me  damping  ratios  (  ~f  ),  on  the  other  hand  vary  erratically 
suggesting  a  non-linear  behavior.  It  can  be  seen  from  equations  (17a), 
(31)  and  (33),  that,  given  u)n  a  constant,  if  -f  is  non-linear,  then  the 
damper  coefficients  are  non-linear. 
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VII.  Conclusions  and  Recommendations 


Conclusions 

Frora  the  analysis  in  Chapter  VI  and  with  the  supplementary  response 
curves  In  Appendix  B,  the-  Vertical  Energy  Absorption  Model  demonstrates 
model  response  characteristics  in  the  roll  and  pitch  modes  that  correlate 
well  to  the  Jindivik  data.  In  the  heave  mode,  results  are  satisfactory, 
but  more  important,  the  potential  for  continued  development  and  more 
precise  prediction  has  been  demonstrated.  The  use  of  a  mechanical 
analog  prediction  scheme  for  an  Air  Cushion  Landing  System  has  sufficient 
merit  to  warrant  further  investigation. 

The  linear  assumption  for  both  the  spring  and  damper  coefficients 
made  in  Chapter  IT  is  born  out  by  analytic  results  and  response  curves 
in  pitch  and  roll.  The  roll  mode  is  superior  in  this  respect. 

Heave  response  of  the  model  is  reasonable,  but  care  must  be 
exercised  in  the  selection  of  the  proper  damping  ratio  (  Y  )  and 
undamped  natural  frequency  (  U),i  ).  For  an  anal;  ical  means  of  finding  Y 
and  ton  ,  the  method  of  averaging  the  values  from  a  peak  time  analysis 
(program  RISE)  provides  good  results.  Reference  the  heave  equation  (8), 
the  terms  on  the  right  hand  side  indicate  a  coupling  of  the  motion  to 
the  pitch  mode.  The  anplicudes  of  the  heave  response  curves  in  Figures 
20a  to  20d  do  not  natch  the  data  consistently  as  -f  and  cOn.  are  varied. 
Using  the  first  peak  values  for  Y  and  results  in  a  good  renuonse 
match  at-  the  first  data  point  (Figure  20a).  but  the  excessive  damping 
ratio  suprcs6es  the  amplitude  of  the  response  at  subsequent  data  poirts. 
Using  the  average  values  of  and  cJo  »  the  response  approaches  thu 
subsequent  data  points  out  overshoots  the  first  data  point.  Those 
observations  suggest  that  the  heave  response  is  non-linear-  The  l-cu- 
linearity  is  in  the  exponential  decay  curve  of  Figure  14  ur.T  if.  d\.  to 
the  product  •  However,  from  the  summary  of  heave  cnar- -torlslics  in 
Table  V,  it  is  shown  that  u)n  is  approximately  constant  and  is  approximately 
equal  to  the  average  value  of  the  heave  With  the  conclusion  that 

the  heave  is  a  com-Luiit ,  tin.-ii  it  lt>  possible  Lu  euneTuue  LV»».  the 
dynamic  heave  spring  coefficient  (-fee  )  is  constant  (linear)  also.  A 
further  observation  is  possible.  With  a  constant  heave  o:,lf  if  the 
product  fuMis  non-linear,  then  the  non-linearity  must  be  due  to  the 
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damping  ratio  ( ~f  ).  Translating  this  observation  to  the  model,  it  is 
now  possible  Lo  conclude  that  the  non-linearity  in  heave  is  due  to  the 
heave  damper  coefficient  (-C-e). 

Of  concern  is  the  fact  that  the  theoritical  response  peaks  of  the 
remaining  pitch  response  after  heave  had  damped  out  did  not  match  the 
data  in  time  (phase  shift)  even  though  there  was  correlation  in  period 
and  amplitude  of  the  exponential  decay.  Since  the  heave  mode  has 
violated  the  original  model  assumption  of  linearity,  it  is  not  possible 
to  form  a  conclusion  about  the  phase  shift  of  the  remaining  pitch  response 
without  tirst  investigating  the  heave  non-linearity. 

In  Chapter  I,  it  was  noted  that  there  are  two  categories  of  trunk 
material,  stretchable  and  non-stretchable.  Since  this  model  investigation 
was  based  on  data  from  an  ACRS  made  of  non-stretchable  material,  use  of 
the  model  with  stretchable  material  trumts  must  be  accomplished  with 
caution.  The  exact  effect  of  a  'tretrl'oble  trunk  on  the  linearity  of  the 
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actcrist^cs  is  net  luncvr.  at  thas  timej  however,  the  model  is 


expected  to  operate  correctly  if  a  proper  expression  for  the  spring  and 
damper  coefficients  are  provided. 

The  schematic  of  the  7EAM  in  Figure  8  visually  suggests  thau  the 
peripheral  springs  rerresent  sections  of  the  trunk  and  that  the  center 
spring  represents  the  cushion.  This  is  not  the  case.  The  model 
requires  a  center  spring  (position  5)  to  account  for  the  heave  dynamics 
in  addition  to  the  roll  and  pitch  linear  springs  which  are  derived  from 
a  torsional  mode  spring.  This  is  illustrated  in  Table  IV  where  the 
linear  s/d  unit  coefficients  are  summarized  from  program  MODANYL* 


Recommendations 

The  following  recommendations  are  derived  from  this  studyi 
1»  An  investigation  into  the  non-linear  nature  of  the  heave 

damping  phenomenon  should  be  made.  A  non-linear  heave  damper 
coefficient  should  be  formulated  and  incorporated  into  the 
model  to  eliminate  the  deficiency  of  the  linear  assumption. 

These  results  should  then  be  used  to  investigate  the  difference 
in  phase  shift  between  the  theoretical  curve  and  the  experimental 
data  of  the  heave  response  after  heave  damps  out. 
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2.  A  method  should  be  investigated  which  would  allow  the  prediction 
of. the  spring  and  damper  coefficients  directly  from  trunK  and 
cushion  parameters  so  that  this  model  can  be  used  to  investigate 
the  dynamics  of  arbitrary  trunk  configurations  without  first 
building  and  testing  on  actual  vehicle  or  scaled  model. 
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Appendix  A 

Analysis  of  JindiviK  Data  (Computer  Solution) 
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Program  DATANYL 

The  prpgrara  for  evaluating  the  Jindivik  ACLS  roll  and  pitch  data 
according  to  the  theory  set  forth  in  Chapter  I V  is  called  DATANYL. 

Table  VI  describes  the  method  for  entering  data  into  the  program  and 
defines  the  symbols  used.  Table  VII  describes  the  output  symbology  and 
units.  Following  Table  VII  is  a  listing  of  DATANYL,  a  single  example 
output  from  the  line  printer,  and  a  plot  of  the  test  data  of  each  of  the 
roll  and  pitch  tests  of  Table  I. 


TABLE  VI 

PROGRAM  DATANYL-INPUT  DATA 

Data 

Group 

No.  of  Data 
Cards  in 
Group 

Format 

Input  Symbol 

1 

n 

X. 

13 

r 

2 

1 

8F5.F10.3, 

A9 

N,P,Q,R,S,W,U, NOTR, MASSEQ, TITLE 

3 

l-N 

FIG. 3 

T 

6 

1-N 

no. 3 

H 

DEFINITIONS  OF  INPUT  SYMBOLS 

F 

N 

P 

Q 

R 

S 

w 

V 

NOTR 

MASSEQ 

TITLE 

T 

H 

The  number  of  sets  of  data  (data  groups  2,3,  £.  4)* 

Total  number  of  data  points 

The  data  point  number  for  the  first  positive  peak  (1  or  2) 
Total  number  of  positive  peaks 

The  data  point  number  of  the  last  positive  peak 

The  data  point  number  of  the  first  negative  peak  (1  or  2) 

Total  number  of  negative  peaks 

The  data  point  number  of  the  last  negative  peak 

The  integer  number  of  the  Jindivik  test  run 

The  moment  of  inertia  (slugs-ft^) 

Enter  "ROLL"  or  "PITCH" 

The  tine  from  t-ofor  all  peaks  in  consecutive  order 

The  adjusted  linear  transducer  heights 

*  data  groups  2,3,  &  4  comprise  a  set  of  data  for  each  Jindivik  test 
run.  Enter  as  many  sets  (F)  ns  there  are  test  runs  to  be  analyzed 
on  a  single  computer  run. 

TABLE  VH 


PROGRAM  DATANYL-OUTPUT  DATA 


Column 

Number 

Column 

Title 

Description 

1 

DATA 

POINT 

The  data  point  number  for  recorded  positive  and 
negative  peaks  in  consecutive  order 

2 

TIME 

(SEC) 

The  time  occurrence  of  the  positive  and  negative 
peaks 

3 

H 

(IN) 

The  recorded  heights  or  perturbations  from 
linear  transducers  for  CG,  Nose,  and  Wing 

4 

T 

(SEC) 

The  period,  measured  between  adjacent  positive 
peaks 

5 

FD 

(CPS) 

The  damped  natural  frequency  in  cycles /second 

6 

WD 

(RPS) 

The  damped  natural  frequency  in  radians/second 

7 

HEQ 

(IN) 

The  equilibrium  height  of  the  system  in  steady 
state 

8 

POSAR 

The  amplitude  ratios  of  the  positive  peaks 

9 

NEGAR 

The  amplitude  ratios  of  the  negative  peaks 

The  following  values 

ace  averaged  to  obtain  the  final  value  used  in  the 

remainder 

of  the  calculations!  T,  FD,  WD,  HEQ,  FOSAR,  NEGAR. 
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Progr.ua  RISK 

Program  RISE  computes  the  namping  ration  (  f  )  and  undamped  natural 
•requency  C  tJn  )  of  the  Jindivik  heave  mode  using  the  method  of  percent 
overshoot  and  rise  time  to  the  first  peak  (Chapter  IV)  and  accounts  for 
an  initial  condition  on  velocity  as  well  as  displacement  at  model 
activation*  The  free  fall  equations  are  included  to  provide  values 
required  by  the  analysis.  To  check  linearity  in  heave,  RISE  can  be 
modified  (lines  381-384)  to  allow  computation  of  •?  and  at  peak  times 
subsequent  to  the  first  peak.  The  program  it  designed  for  use  on  the 
teletype  remote  terminal.  Table  VIII  provides  the  method  for  reading 
data  into  the  program.  Table  IX  describes  the  symbols  and  units  for 
output.  Symbols  in  both  Tables  are  in  accord  with  Figure  IS.  Figure  31 
presents  a  listing  of  the  program  RISE.  A  summary  of  all  calculations 
from  program  RISE  can  be  found  in  Tables  II  and  III  of  Chaper  IV. 


table  VIII 

PROGRAM  RISE-INPUT  DATA 

Input 

Statement 

Format 

Input  Symbol 

Read  1 

5F10.3 

TAU  TAUPI  H I  HEQ  HPI 

Read  4 

2F10»4,2I5 

ZETA  WN  N  L 

DEFINITIONS  OF  INPUT  SYMBOLS* 

TAU 

TAUP1 

HI 

HEQ 

HP1 

ZETA 

WN 

H 

L 

Time  of  free  fall  (sec.) 

Time  from  release  to  the  first  peak  (sec.) 

Height  of  release  point  from  surface  (in.) 

Equilibrium  height  from  surface  (in.) 

Height  of  first  peak  from  surface  (in.) 

Damping  ratio  -  assumed  or  guess  for  iteration  (N/D) 

Undamped  natural  frequency  -  assumed  or  guess  for 
iteration  (rad. /see.) 

Integer  value  of  cycles  -  for  fii.su  peak 

Integer  value  of  number  of  iterations  desired 

*  Symbols  arc  in  accordance  with  Figure  16,  Chapter  IV 


TABLE  IX 

■ 

PROGRAM  RISE  -  OUTPUT  DATA 

Output 

Statement 

Output 

Symbol 

Output  Definition 

Units 

Print  11 

TP 

Time  to  first  peak 

sec* 

20 

Initial  displacement  from  equili¬ 
brium,  positive  measured  down 

in. 

Z0D0T 

Initial  velocity  at  to 

in. /sec. 

2  PI 

Displacement  of  the  first  peak 
from  equilibrium 

in. 

Print  15 

N 

The  N  value  read  in 

N/D 

ZETA1 

The  first  guess  for  ZETA 

N/D 

WN1 

The  first  guess  for  WN 

rad. /sec. 

Print  13 

ZETA 

The  iterative  solution  for  ZETA 

N/D 

WN 

The  iterative  solution  for  WN 

rad. /sec. 

62 
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Pro pram  HODANYL 

The  evaluation  of  the  equations  of  notion  of  the  Vertical  Energy 
Absorption  Model  (VEAM)  in  its  three  degrees  of  freedom  (DOF)  is 
Accomplished  by  program  MODANYL.  A  listing  of  the  program  is  presented 
In  Figure  32.  Due  to  program  volume,  only  the  graphical  results  are 
included  in  this  study. 

Thie  program  is  designed  to  solve  the  VEAM  equations  with  the  use  of 
AFIT  Subroutine  RKDES ,  a  differential  equation  solver  using  r  v  riable 
•tep  fourth-order  Runga-Kutta  method.  In  addition,  the  mode  ; -rings 
and  dampers  arc  converted  to  linear  spring/daraper  units  with  the 
equations  developed  in  Chapter  V- 

The  total  vertical  displacement  as  a  result  of  movement  in  the  3 
DOF  of  any  arbitrary  point  in  the  plane  of  the  model  is  computed  by 
Subroutine  LOCATE.  Displacement  of  any  point  (P)  relative  to  the  CG 

could  be  accomplished  by 
LOCATE,  but  the  model  assumptions 

Cf-  r  **  TT  »  «  .  ■  "■  i  W  1  I  m!  fit*. 

*  *•  -a  ^  a. t  AA<uAWAil^ 

the  arbitrary  point  to  the 
horizontal  plane  (Z*  »  0). 

For  this  study,  the  three 
positions  in  the  xy  plane  that 
coincided  with  the  locations  of 
the  Jindivik  linear  displacement 
transducers  were  chosen  for 
output  in  order  to  compare  with 
the  test  data. 

A  subroutine  called  ANAILT 
plots  the  theoretical  curve  and 
the  experimental  data  points  for 
comparison.  Since  this  plot  routine  is  similar  to  DATAPLT  of  Appendix  A, 
a  listing  is  omitted. 


Small 


Fig.  32.  Vertical  Displacement 
for  Subroutine  LOCATE 
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Th«  equations  of  motion  solve  a  dynamic  force  balance  relative  to  the 
equilibriura  axis.  All  initial  conditions  on  displacement  that  are  applied 
to  the  equations  of  motion  must  be  referenced  to  the  equilibriura  axis  also* 
Collecting  data  relative  to  the  equilibrium  axis  is  difficult  aince  the 
equilibriura  position  changes  according  to  test  configuration.  Therefore, 
data  is  measured  from  a  known  horizontal  such  as  the  surface  or  the 
spring-neutral/contact  plane.  For  the  purposes  of  these  calculations,  the 
spring-neutral  or  contact  plane  is  chosen  as  the  reference.  A  total 
force  and  moment  balance  is  accomplished  and  the  equilibrium  displacement 
of  the  model  from  the  spring  neutral  plane  is  calculated  for  each  mode. 

The  following  relationships  provide  the  means  to  reference  the  initial 
conditions  to  the  equilibrium  axis. 


W  ^ 


£  W 


•*c  Z  Si. ,  tZ/S1 

'<■»  V  *./ 


<t>  *  J«~  h 

-  Vi 


0  for  all  timesj  that  is, 


05) 

(36) 

(37) 


the  equilibrium  axis  and  the  contact  axis  measure  roll  relative  ro. 
the  horizontal. 
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Fig.  34.  Pro&ran  MODANYL-,  Subroutines  F  and  LOCATE 


EXPERIMENTAL 


EXPERlflENTRL 


Fig.  37a.  Model  Response,  Heave  Teat  121,  CG  Transducer 


<$ 


Nose  3°  below  equilibrium 


EXPERIMENTAL 


Fig.  39.  Model  Response,  Pitch  Test  59,  Nose  Transducer 
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Appendix  C 

Justification  of  tho  Snail  Anp,le  and 
Horizontal  Displacement  Assumptions 
in  Pitch  and  Roll 
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General 

In  tl. -  ma in  text,  simplifications  were  made  to  the  pitch  and  roil 
analysis  that  were  based  upon  the  assumption  that  the  resect ive  angles 
were  small  and/or  the  assumption  t>at  the  horizontal  displacement  of  a 
point  during  pitch  or  roll  was  negligible.  Justification  of  these 
assumptions  in  the  text  would  have  been  unwieldy,  ajid  would  have  detracted 
from  the  primary  analysis.  The  roll  mode  has  been  chosen  for  thi3 
detailed  illustration  because  it  has  the  maximum  motion. 

Geometry 

Tne  roll  linear  displacement  transducer  was  affixed  to  the  Jindivil; 
suspension  frame  a  distance  CL)  above  the  horizontal  surface  through 
the  center  of  gravity.  A  vire  on  inertia-reel  from  the  transducer 
was  attached  to  the  wing  tip  at  point  a.  Figure  41  shows  the  wing  in 
deflection  from  equilibrium  (  ).  The  transducer  record^  a  vertical 

displacement  (A),  but  the  attachment  point  is  constrained  by  the  wing 
and  must  sweep  out  the  arc  w-w.  The  wire,  oric.inaliy  °f  length  (  L  ),  is 
reduced  to  a  length  (L*A)«  Point  a,  attached  to  the  wire,  muse  be  or 
arc  t-t  also,  where  arc  t-t  is  jxirt  of  a  circle  centered 
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at  the  transducer  with  radius  (L-Ji;).  If  the  movement  were  truly 
vertical,  the  point  of  attachment  would  be  at  a*,  but  it  is  at  position 
&  instead.  This  results  in  a  vertical  error  (A.^t),  and  a  horizontal 
error  (  Xy^  -  X<x  )  • 


Analysis 

An  origin  is  created  at  the  center  of  gravity  with  an  x,y  axis 
defined  for  this  problem  as  illustrated.  Equations  are  then  written 
which  describe  the  arc3. 

For  arc  w-wi  ?.1  -t  31-  =.  (3S) 

For  arc  t-t.  “  W*  *  (  '  i-Y’  =  (L"^  (39) 

Simultaneous  solution  of  these  equations  will  yield  the  coordinates  for 
points  a  and  f.  From  (3S),  solve  for  X  =  ^  and  substitute 

into  (39)  to  solve  for  _ _ _ _ — . 

LK  -  luM*  -  M [  L 1  -  Pv/jL  (\)''  ) 

= - - - -  <40> 


where 


a  [  l*  .(^y-  j 


K  =  2(l>/i)1.2AL-A‘-  (41) 

r 

For  the  maximum  roil  case  y 

L  *=  14  ft.  “  168  in. 
tyz  «=  164.75  in. 

A  »  22  in. 


from  which 

ba=  22.00b  in. 

.  n  _ _ 

,uub  in. 

This  calculation  shows  that  the  vertical  error  is  negligible. 
Substituting  into  equation  (38)  results  in  a  solution  for 
X**  163.268  in. 

Point  a  translates  by  the  amount  Vz  ~  Xa  c  1.482  in.  Thus ,  the 
horizontal  translation  is  only  about  11  of  the  half  span  (  b/2 -)»  and  can 
thus  be  considered  negligible. 

The  method  above  allows  a  precise  calculation  of  the  roll  angle, 
also,  where 

*  .  .  1  u„  / 

(j>  -  too.  /  Xo, 

For  roll,  -  7.08  degrees,  which  supports  the  small  angle 

assumption. 

The  pitch  case  follows  the  same  analysis  except  the  half-span  is 
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replaced  by  the  distance  between  the  nose  transducer  attachment  point 
and  the  center  cf  gravity  (166.51  in.),  and  the  m^inum  vertical 
displacement  is  19.80  inches.  The  analysis  shows  that 
,ka  19.8047  in. 

165,328  in. 

&5ih  .00n7  in. 

G*«=  6.83  degrees 
Translati.cn  =  1.182  in. 

vhich  together  support  the  original  assumptions. 

Although  the  above  analysis  used  the  displacement  transducers  as 
the  example,  the  same  assumptions  hold  true  for  the  trunk.  The  x- 
displacement  changes  magnitude  but  their  relationships  and  the  maximum 
angles  remain  the  same.  The  results  of  this  analysis  were  used  during 
development  of  the  ACLS  model  as  described  in  Chapter  II  and  in 
subsequent  calculations. 
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Appendix  l) 

AFFDL  Jit.iUvik  and  Test  Data 


ga;i/al/,_v-i 


Certain  data  required  in  various  conputat iona  was  provided  by 

"v 

v  }  the  Air  Force  Flight  Dynamics  Laboratory,  Mechanical  Dranch. 

1.  Ejector  Tost  Results 


Equilibrium  Hover 

Conditions 

Vent 

Opened 

Vent 

Closed 

Flow  Added 
P<L  -12  ?psi 

Flow  Added 
Pc-3c.4  Psi 

Trunk  pressure  (  Pt  ) 

1.94 

1.70 

1.64 

1.56 

Cushion  pressure  C  Pc  ) 

0.U3 

0.61 

0.65 

0.68 

Trunk  flew  (mt)  lbm/sec 

0.79 

0.97 

1.02 

1.07 

Cushion  flow  (eu.)  lb^sec 

UNTC 

0.0 

0.46 

1.32 

Trunk  footprint  (  Pt  ) 

00 

n 

(V 

r  * 

617 

579 

565 

2.  Static-  Stiffness 


Roll  (JJlft)  ft~lb/rad 

45S3.66 

2437.75 

2864.79 

- 

Pitch  (itp)  ft-lb/rad 

51,566.2 

51,566.2 

51,566.2 

- 

Heave  (Ah  )  lb/ft 

9000 

8400 

7  203 

- 

Heave  Test 

Drop  Tine, 

Drtl  «*  n  f-r\  4-  /-.vn  f*  n  r»  •*  ( 

0.1 

” 

0.093 

0.081  ] 
I 

Jindivik  and  ACRS 

Moment  of  Inertia  (Jxx) 

slug-ft^ 

1190 

Moment  of  Inertia  (Jyy) 

,  ..  "> 

S lUg-tt 

1810 

Weight 

lbf 

2470 

Trunk  width 

in 

35.2 

Trunk  length 

in 

107.7 

CP  to  CG  Displacement  (S) 

in 

8.56 
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